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Background: Models to calculate small isospin-symmetry-breaking effects in superallowed Fermi decays have 
been placed under scrutiny in recent years. A stringent test of these models is to measure transitions for which 
the correction is predicted to be large. The decay of '^^Cl decay provides such a test case. 

Purpose: To improve the 7 yields following the /? decay of "^^Cl and to determine the ft values of the the fi 
branches, particularly the one to the isobaric-analogue state in ^^S. 

Method: Reaction-produced and recoil-spectrometer-separated ^^Cl is collected in tape and transported to a 
counting location where /3 — 7 coincidences are measured with a precisely-calibrated HPGe detector. 

Results: The precision on the 7 yields for most of the known /3 branches has been improved by about an order of 
magnitude, and many new transitions have been observed. We have determined '^^Cl-decay transition strengths 
extending up to ~ 11 MeV. The ft value for the decay to the isobaric-analogue state in ^^S has been measured. 
A comparison to a shell- model calculation shows good agreement. 

Conclusions: We have experimentally determined the isospin-symmetry-breaking correction to the superallowed 
transition of this decay to be {5c — 5Ns)exp = 5.4(9)%, significantly larger than for any other known superallowed 
Fermi transition. This correction agrees with a shell-model calculation, which yields 5c — iSns = 4.8(5)%. Our 
results also provide a way to improve the measured ft values for the /3 decay of '^^Ar. 

PACS numbers: 23.40.Bw, 24.80.+y, 29.30.Kv, 23.20.Lv 



MOTIVATION 

The comparative half-lives of superallowed Fermi /3 
decays between 0+ isobaric analogue states have been 
the focus of intense research activity for many years and 
presently represents one of the most stringent tests of 
the Standard Model of the electroweak interaction [ij. 
The high precision of both experimental measurements 
and theoretical calculations of their ft values set strin- 
gent limits on scalar and right-handed currents, verify 
conservation of the vector current to ~ 10~*, and de- 
termine the up-down element of the Cabibbo-Kobayashi- 
Maskawa (CKM) quark-mixing matrix, Vud [IHll ■ Exper- 
imentally, the ft value of thirteen cases have been mea- 
sured to < 0.3%; this places a demanding requirement on 
theory to attain similar precision. Although these tran- 
sitions are intrinsically simpler to describe theoretically 
than most (3 decays because they are relatively insensi- 
tive to nuclear-structure effects, small (~ 1%) corrections 
must be applied to account for the fact that the decay 
occurs within the nuclear medium. Recently, emphasis 
has been placed on scrutinizing the nuclear-structure- 
dependent isospin-symmetry-breaking (ISB) corrections, 
5c 041] I which characterizes the degree to which the 
Fermi matrix element, Mp, deviates from Aig, its value 
in the limit of strict isospin symmetry: 



The 13 most precisely-measured cases mentioned pre- 
viously are all isospin T = 1 to T = 1 transitions in 
A = 4n + 2 nuclei. Shell-model calculations for these 
cases yield values of order Sc 0.5% for ^ < 56 [l| and 
values of order Sc - 1.5% for ^ > 56 [l|, l^. If attention 
is switched to ^ = 4n nuclei, even larger values of Sc 
are predicted, which if experimentally extracted, would 
provide an even more demanding test of such ISB calcu- 
lations. The reason larger values are expected in A = An 
nuclei is that the daughter analog state sits among many 
states of lower isospin, T — 1. Some of these states have 
the same spin as the analog state and sizable isospin mix- 
ing can occur. 

In this work, we will expand on a recent Letter [lol | 
which discusses an extraction of the isospin-symmetry 
breaking correction in the decay of 1+,T = 1 '^^Cl. Its 
Fermi decay branch feeds a 1+, r = 1 state in '^^S, whose 
position in the spectrum at 7001-keV excitation [11| is 
very close to a known 1+,T = state at 7190 keV |12j . 
As discussed below, a calculation of 6c for this case yields 
Sc = 4.8(5)%, a significantly larger value than those 
found in A = An -\- 2 nuclei. 

Another motivation for this work is related to the re- 
cently measured T = 2 decay of '^^Ar [l3|. Here calcu- 



lations predict Sc = 1.8% [14|. The ISB correction for 



Wf^" = \Mo\''{l-Sc). 
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this case extracted from the experimental ft values was 
(2.1 ± 0.8)% [11 and later corrected 
(1.8 ± 0.8)% in Ref. [15|], where an improved 
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value for the end-point energy was deduced. A poten- 
tially large source of systematic uncertainty arising from 
the need to detect 7 rays in this measurement may be 
minimized using the 7 branches of ^^Cl. This is because (3 
decay of ^'^Ai is followed by the /3 decay of ^201 64.4(2)% 
of the time Thus, ^^Cl provides an in situ efficiency 
calibration which is useful to extract isospin-breaking in- 
formation from ■^^Ar. The present work opens the possi- 
bility for significant improvements in the precision with 
which 5c can be determined in the decay of '^^Ar. 



EXPERIMENTAL PROCEDURE 

The experiment was carried out at the Cyclotron In- 
stitute, Texas A&M University. A primary beam of ^'^S 
was produced by an ECR ion source and injected into 
the K500 superconducting cyclotron to accelerate it to 
K. 24.8 MeV/nucleon. The « 400 nA ^'^S beam ex- 
ited the cyclotron and was directed towards the target 
chamber of the Momentum Achromatic Recoil Separa- 
tor (MARS) [Hi. A 20 MeV/nucleon secondary beam 
of CI was produced via the inverse kinematic transfer 
reaction, ^Hi^'^'&.nY'^Ci on a LN2 cooled, hydrogen gas 
target at « 1.4 atm. MARS was used to spatially sepa- 
rate the reaction products, resulting in a ^^Cl beam with 
an intensity of 2 x 10^ ions/s. Beam contamination was 
identified using a position-sensitive Si-strip (A_E) detec- 
tor followed by a silicon {E) detector which were placed 
just downstream of the MARS focal plane. For data col- 
lection, the Si detectors were removed and the beam ex- 
ited the MARS beamline through a 50 ^m-thick Kapton 
window. The beam then passed through a 0.3 mm thick 
BC404 scintillator to count the number of ions. Prior 
to being implanted into a 76 ^m-thick aluminized-Mylar 
tape which is part of a fast tape-transport system, the 
beam was passed through a set of Al degraders. The 
thickness of the degraders was chosen to ensure that the 
activity was deposited mid-way through the tape. The 
different ranges of the contaminants compared to ^^Cl al- 
lowed further purification; however, since we were search- 
ing for small branches and wanted to maximize the yields, 
we allowed greater contamination than usual of the de- 
posited activity, accepting 91% as our final purity of '^^Cl. 

The ^^G\ atoms were collected in an « 1 cm diame- 
ter spot on the tape for 0.8 s, after which the beam was 
interrupted and the tape-transport system was triggered 
to move the activity to a shielded counting station 90 cm 
away. The latter was accomplished in « 180 ms. The set- 
up is shown schematically in Fig. [TJ Once transported to 
the shielded area, /? — 7 coincident data were acquired 
for typically 1 sec (83% of the total data set). In a few 
of the runs (corresponding to 11% and 6% of the data 
respectively) we used count times of 2 sees and 4 sees to 
check for long-lived contaminants. The data were regis- 
tered event-by-event by recording all /3 — 7 coincidences 



MARS 
beamline 




aluminized 
Mylar tape 



implantation 
tape-transport 
count decays 



< ] 


L c 


/cle >| 





























0.80 0.98 1.98 2.78 2.96 

t [sec] 

FIG. 1. Schematic diagram of the end of MARS and the fast 
tape-transport system as well as a timing diagram for the 
experiment. The counting time for decays was generally 1 sec 
as indicated, but times of 2 and 4 sec were also used for the 
diagnosis of longer-hved contaminants. 



between a 1-mm-thick BC404 plastic scintillator and a 
70% HPGe detector. The 1.5 inch diameter scintillator 
AE detector had a threshold of 40 keV and was placed 
5 mm behind the tape subtending « 32% of the total 
solid angle for the /3s. The 7 detector was placed much 
farther away: 15.1 cm from the tape to reduce the effects 
of coincidence summing. The 7-ray energy, the /S.E of 
the /3, the coincidence time between them, and the time 
of the event relative to the beginning of the cycle were all 
recorded. The typical tape cycle, outlined in Fig. [2 was 
repeated continuously throughout the experiment. The 
total number of /3 singles events and the total number of 
heavy ions (His) from MARS (detected by the first scin- 
tillator) for each cycle were determined from scalers and 
recorded. The ratio of /3 singles to His was used to veto 
bad cycles where, for example, the tape transport did not 
place the activity exactly in the correct location. Of the 
« 36 000 cycles made over the course of the experiment, 
92.5% survived the /3/HI ratio cut. 



7 Efficiencies 

The extremely precise absolute photopeak efficiency 
calibration of the HPGe detector is a critical compo- 
nent of this experiment. This efficiency has been care- 
fully studied up to energies of 3.5 MeV as discussed 
in Rcfs. [i3-[li- Over the range of 50 - 1400 keV, 
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FIG. 2. Top: the adopted absolute photopeak efficiency curve 
of the HPGe used in this worfc. The solid line is the mea- 
sured efficiency up to 3.5 MeV from Refs. [1343 > while the 
dashed line above 3.5 MeV represents an extrapolation based 
on Cyltran calculations. Bottom: percent difference of a 
Penelope simulation minus the adopted curve (dotted line) 
and the assigned uncertainties to the efficiency curve used in 
the present analysis (hatched region). 
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FIG. 3. Total efficiency of the HPGe as simulated by Pene- 
lope. The filled circles represent the simulated points with 
our basic geometry (detectors and Mylar tape) while the open 
circles show the results when our basic geometry is encased in 
an Aluminum cylinder. The dashed lines going through each 
is a 4**^ order polynomial fit, and the solid central line is the 
average of the two which we take to be our total efficiency 
curve. We assign an uncertainty that spans the result of both 
simulations. We used these curves to determine systematic 
uncertainties associated with summing in our detectors. 



Rcf. [l8| discusses how the absolute efSciency has been 
calibrated to ±0.2%, and Ref. [l^ uses a combination of 
measurements and calculations using the Cyltran [i^l 
Monte Carlo code to extend the efficiency curve from 
1.4 - 3.5 MeV to ±0.4% precision. For this work, the 
highest energy 7 rays observed are at 7.2 MeV, requiring 
us to further extend the photopeak efficiency curve. Our 
adopted curve and its extrapolation, shown in the top 
panel of Fig. [51 is calculated using the the same Cyltran 
program used in Refs. 17-1^. Note that above 3.5 MeV, 
the efficiency curve is no longer linear on the log-log plot 
as a simple extrapolation might suggest, instead falling 
off more quickly. We additionally checked this extrapo- 
lation against an independent calculation based on the 



Monte Carlo radiation transport code Penelope [21 



As can be seen in the bottom panel of Fig. [2l there is 
excellent agreement between the two calculated photo- 
peak efficiencies over the range of energies observed in 
this work. As the figure also indicates, we increase the 
uncertainties in the efficiency curve, adopting conserva- 
tive uncertainties of ±1% from 3.5-5 MeV, and ±5% 
above 5 MeV. The differences between the Cyltran and 
Penelope extrapolated efficiency curves are well within 
these uncertainty ranges. 

We have investigated the effects of summing in the 
HPGe detector using our Monte Carlo simulations. A 
small but non- negligible factor arises from Compton sum- 
ming of 7 rays that scatter off various volumes. This 
makes knowing the total efficiency for 7 detection nec- 
essary, although the precision does not need to be very 



stringent because this summing with photopeak events 
is a small correction. However, accurately quantifying 
the effects of summing is difficult due to the large vol- 
ume needed for tracking photons that may scatter from 
any of the surrounding elements. Our simulations where 
the geometry contained only detailed descriptions of the 
detectors and Mylar tape (i.e. neglecting the table, the 
floor, the walls, etc.) underestimates this summing ef- 
fect. Rather than attempting to include the many po- 
tentially important scattering surfaces, we chose to per- 
form Penelope simulations with our geometry encased 
by an Aluminum cylinder knowing that this surely over- 
estimates the effect. The result of simulations using these 
two geometries is shown in Fig. [3] where differences as 
large as a factor of two arise at lower energies. We take 
the total efficiency to be halfway between these two simu- 
lated curves, with an uncertainty that spanned the results 
of both. This results in a large uncertainty in the total 
efficiency but since the Compton summing with photo- 
peaks is a small correction, this conservative estimate 
does not limit our determination of the 7 branches. 

We also investigated the possible effects of 7-ray an- 
gular correlations, which may affect the probability of a 
photopeak event summing with the photopeak from an- 
other cascade 7. The only cases that would result in a 
non-zero 7-7 correlation are 1+— > 2+— > 0+ cascades; we 
tested both E2 and Ml transition types and found that 
independent of the transition type, 7-7 angular corre- 
lations do not lead to any additional summing in our 
geometry. 
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Qec and the /3 Efficiency 

The mass excess of '^^Cl is obtained from our averaging 
Refs. [mil,!!! to get ME(32C1) = -13334.60(57) keV, 
and ME(32S) = -26015.535(2) keV is taken from 
Ref. [3]. Taken together, the decay energy is Qec = 
12680.9(6) keV. Although the fraction of events below 
the finite threshold of 40 keV in the plastic scintil- 
lator is expected to have a weak dependence on the 13 
end-point energy, in principle the /3 efficiency depends 
on Q/3- To investigate this effect, Penelope simulations 
were used to determine the /3 efficiency of the plastic 
scintillator. A linear fit of the mean /3 efficiency for the 
entire (3 spectrum versus /3 end-point energy yielded an 
efficiency that was consistent with being flat over the 
4 — 12-MeV range of end-points relevant for this decay: 
rj{Qp) = [32.32(41) - 0.022(45)Q/3]%. 



DATA ANALYSIS 

Figure 13] shows the spectrum from the resistive-readout 
position-sensitive Silicon detector at the focal plane of 
MARS. This logarithmic 2D plot shows that the most 
significant contaminations in the beam prior to our clos- 
ing the purifying slits were ^°S and "^^S, with ^^Cl making 
up « 86% of the beam at the focal plane of MARS. The 
^^S contamination was minimized by closing the vertical 
slits at the focal plane of MARS as indicated in the fig- 
ure, reducing it from « 3% to « 0.4%. The slits had 
little effect at removing the '^'^S contamination, however, 
as it lies in the same vertical band as the '^^Cl. The pu- 
rity of the beam at the focal plane of MARS with the 
vertical slits in place was ~ 89%. As mentioned previ- 
ously, the 30-31 S contaminations were further reduced by 
the degraders, which were chosen to maximize the im- 
plantation of ■^^Cl in the centre of the Mylar tape. The 
purity of the beam implanted in the tape was improved 
by another couple of percent to w 91% ^^Cl with ^°S the 
largest contaminant at « 7.5%. Note that we could have 
obtained a higher purity at the cost of reducing the rate; 
however, since the 7 energies of the contaminants do not 
overlap the ^^Cl lines, we chose to maximize the rate. 

The plot in Fig. [5] shows the time difference between /3 
particles detected in the scintillator and 7-rays observed 
in the HPGe detector plotted against the 7 energy. One 
clearly sees a strong peak at tp—t^ « 475 ns. In addition 
to the expected walk at low 7 energies (below 511 keV), 
a smaller 2"^^ peak around 600 ns in the timing is visi- 
ble; although the source of this peak is not fully under- 
stood, the 7 spectrum gated on just it proves that these 
are good events. We therefore included it in the analy- 
sis and defined a time gate between 430 and 800 ns to 
select real coincidences (the dashed horizontal lines in 
Fig. [5]). The rest of the timing spectrum defines another 
window which selects accidental coincidences. The Er^ 
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FIG. 4. Identification of the ion beam at the focal plane 
of MARS, before the purifying slits in the vertical position 
were moved to the closed positions (indicated by the hatched 
areas). The region labeled ^^Cl represents 86% of all events 
in the spectrum, which was increased to 89% once the shts 
were in place. The tails in the vertical position are a result 
of incomplete charge collection in the resistive-readout Silicon 
/S.E detector. 



projection of Fig. [S] shows two curves: the top one repre- 
sents the projection of the real-coincidence window, and 
the lower one represents the accidental coincidences. The 
lower one, properly normalized, was then used as a back- 
ground spectrum to be included in our fitting function as 
described below. The real-coincidence window position 
and size were varied significantly to check for potential 
systematic errors. This procedure yielded results with 
no significant sensitivity to the particular window used, 
as long as it covered the range containing all of the real 
coincidences. 

Figure [6] shows the 7 spectrum observed in the HPGe 
detector in hardware coincidence with a (3 signal in the 
scintillator. This is the same as the upper projection of 
Fig.[5]but with finer binning. Except for a strong peak at 
677 keV from the decay of "^^S nearly every statistically 
significant peak is associated with the decay of ^^Cl. One 
exception is at 2776 keV where 360 ± 50 counts are seen 
which could not be identified based on the known levels 
in ^^S 12[ nor with any contaminants. If this was, in fact, 
a product of the decay of '^^Cl, it would only represent a 



0.1% 7-ray yield. 

A reassuring check of the cleanliness of our data and 
identification of the '^"S contamination is shown in Fig.[7l 
This shows a comparison of the lifetime of the two most 
intense peaks in the '^^Cl spectrum (2230 and 4771 keV) 
as well as the main peak (677 keV) from ■^°S. Dead time 
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FIG. 5. Logarithmic 2D plot of the 7 energy versus the 
timing between the /3 and 7. Accidental coincident events 
were defined to have tp-y < 430 ns or tp-/ > 800 ns. The 
lower curve of the projection of the HPGe energy repre- 
sents these background events (scaled according to the cut- 
window size). The upper projection is the good events with 
430 ns < tp-y < 800 ns. 



and pile-up effects were assumed to be negligible cor- 
rections and are not included in these half-life curves. 
The fit was from 0.060 to 0.985 s after the activity was 
transferred to the counting station. The fit lifetime of 
2230-kcV 7 events is 0.3012(13) s which is consistent with 
the 0.300(5) s lifetime fit from 4771-keV events. Both of 
these results are in a gree ment with the accepted half-life 
of 0.298(1) s of Ref. For the 667-keV line from 
events, the Compton tail from higher-energy 7s from the 
shorter- lived ^^Cl represent a contamination to the '^°S 
curve. In order to remove their effect, the minimum time 
included in the fit range for '^°S was 1 sec, over three 
■^■^Cl half-lives. The fit yielded the half-life in this case 
to be 1.169(34) s, in good agreement with the accepted 
value of 1.1786(45) s from Ref. Other ^"^Cl peaks 

did not have enough statistics to be able to confirm the 
consistency of their tif2 decay curves. 

The 7 Peak Areas 

In order to extract peak areas we used a fitting function 
consisting of four terms: 



1. A Gaussian, Foa 



corresponds to a (5-function at the 7 peak chan- 
nel number, /x, convoluted with noise of width a. 
Due to the sharp rise of this peak, we integrated 
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FIG. 6. The 7 spectrum observed by the HPGe detector in 
coincidence with a /3. Peaks are labeled with their energy, 
the symbol * or the symbol **, referring to the full energy, 
single-escape and double-escape peaks, respectively. Peaks 
not associated with ^^Cl have in addition a label indicat- 
ing the parent nucleus. Decays from the 7001-keV isobaric 
analogue state are highlighted with boxed values. The only 
significant background peak is at 677 keV from the '^"S con- 
tamination. 



the function over the ADC channel width. Ax, be- 
fore comparing to the experimentally observed his- 
togram. Thus the main peak of our fitting function 
is: 



-fj.-Ax/2 
V^cr 



Fix) = - 

where erf is the error function. 

2. A convolution of an exponential, Fcxp = 
exp [—X(x — /x)], with a Gaussian. Normalizing to 
unit area over the range zero to infinity, we define 
a low-energy tail function corresponding to incom- 
plete charge collection as: 



Tix) 



A exp + A(a; - /i) ) erfc 



( x—fi+Xo 
V2a 



l + ^rf(^j-^-P(^-V)erfc^^ 
Here erfc(x) — 1 — erf(a;) is the complementary 
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FIG. 7. Number of counts as a function of time after the activ- 
ity was transferred by the fast tape-transport system (points) 
and Ufetime fits (sohd hnes). These decay spectra are gated 
on the two most prominent ^^Cl hnes (2230 and 4771 keV) 
as well as the dominant contaminant, at 677 keV. The 
^^Cl peaks were fit to runs where events were counted for 1 s 
following transfer, and the longer-lived peak used runs 
where the count time was 4 s. 

error function and A is the decay constant of the 
exponcntiaL 

3. A constant background. -Fbkgd = b. 

4. The background histogram from the accidental co- 
incidences in /3 — 7 timing, Faccid(2:), as described 
previously. 

In order that the total number of 7 signal events is a 
free parameter of the fits, we combined and normalized 
the terms so that our final fitting function was: 

Ffit = [(1 - al,,)F{x) + al,,T{x)] 

+ Fbkgd +i^accid(x). (2) 

In our model, the parameters and ataii are squared 
so that neither the normalization of the photopeak area 
nor the relative amplitude of the tail due to incomplete 
light collection are able to converge to a negative value 
when fitting very small amplitude or non-existant peaks. 
The total number of good 7 events in a peak is given by 
N'' = and the uncertainty is AA^''' = 2a^Aa^, where 
Aa^ is the statistical uncertainty returned from the fit- 
ting routine. Note that by defining the fitting function 
in this way, our statistical uncertainty in the number of 
7 events includes any correlations with other parameters 
of the fit. 

The 7 energies investigated included any transitions 
between states such that > 400 keV. The background 
from the Compton edge of the 511 keV peak compromised 
the sensitivity of searching for peaks below this energy. 
All fits were made using a FORTRAN code based on the 
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FIG. 8. (Colour online) The top panel shows a typical fit, in 
this case to the 2230-keV peak, using the fitting function of 
Eq. The dashed line is the tail, T{x), from incomplete 
charge collection; the dotted is the accidental coincidences, 
faccid; the dash-dot line is a flat background; and the solid 
histogram is the total fitting function. The bottom panel 
shows the residual of the fit in terms of a. 

Marquardt algorithm (27| for minimization. The pro- 
gram assumes Poisson statistics in the data (rather than 
Gaussian), and thus properly handles bins with very few 
counts. The data were divided into small blocks with 
equal number of total counts and analyzed individually. 
This procedure minimized the effects of gain variations 
since each data set was acquired over a limited range of 
time rather than over the whole run. Gain variations 
determined by the change in the fit centroid value were 
always below 0.05%; so in the end, the effects due to gain 
variations were negligible. The number of blocks of data 
varied for each peak depending on its intensity; for ex- 
ample, while fitting the high- statistics 2230-keV peak the 
data were divided into 170 blocks, while the much weaker 
7189-keV peak had all of the data summed together into 
one histogram before fitting. When the data were divided 
into blocks for fitting, the areas were summed and un- 
certainties propagated to get the total number of counts. 
We note that this number is in statistical agreement with 
the result obtained when fitting the summed histogram, 
although in the latter case the was worse due to small 
drifts in the gain. 

Figure [8] shows a typical fit to the 2230-kcV peak. In 
general the fits were excellent, with very few converging 
with a confidence level (CL) far from the ideal 50%; the 
distribution of confidence levels has a mean of 30.8% and 
a standard deviation of 24.2%. 

The energy calibration of the energy of our HPGe 
detector was made using: (a) the two most precisely 
known and strongest 7 de-excitations from levels in '^^S 
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(2230.49(15) and 4281.51(26) keV from Ref. and 
(b) the three independent 7-rays from our main contam- 
inant (667.01(3), 708.70(3) and 2342.2(1) keV 



The /3 and 7 Branches 

Another FORTRAN code, also based on the Marquardt 
algorithm, was used to fit /? and 7 branches to the ob- 
served number of 7 events, N^j^ where j represents the 
state to which the z*** state 7 decays. Letting f3i repre- 
sent the probabihty for the (3 decay to proceed to state 
i, we normaUze the B branches such that 



(3) 



where i is summed over all of the states considered in 
the analysis (see below). Similarly, we normalize the 7 
branch probabilities, Jij, from a given '^^S excited state 
i to lower state j as 



1. 



(4) 



In terms of these branches, to first order the number 
of observed 7 rays is given by: 



Nu 



k>i 



(5) 



where the first term in the brackets arises from the /3 
transition directly to state i, and the second term repre- 
sents the feeding from /3 transitions to higher levels that 
7 decay to state i. Here Nt^t is the total number of de- 
cays (a free parameter in the fits since the total number 
of /3's was not precisely measured), e^j- is the photopeak 
efficiency of the HPGe detector at an energy Ei,j, and 
T]i is the efficiency of the plastic scintillator to observe 
a P with an end-point energy Qi. Small corrections to 
Eq. ^ (included in the analysis but omitted here for 
clarity) are required to account for (a) summing with 
cascade 7s from above and below (which requires the to- 
tal efficiency of the HPGe), and (b) summing with the 
511 annihilation radiation since this is a /?"*" decay. 

The nine lowest (0, 1, 2)"*" states observed by Detraz et 
al. ^ and Anderson et al. ^ along with the (1.0l^ j)% 
ground state branch from Armini et al. [25| represent 
most of the /3 yield; however little is known above 
7.2 MeV of excitation energy. Some /3-delayed proton 
and alpha decays have been observed [31| to states above 
8.7 MeV with particle energies as low as 762(5) keV, but 
this still leaves a 1.5 MeV window of Q-value in which 
no /3 transitions have been identified. This suggests that 
there is no appreciable (3 feeding in this energy region; 
however, it does not rule out the possibility of a large 
number of weak /? transitions. Each of these transitions 
may be too small to be detected individually, but could 



cumulatively contribute a total /3 strength of up to a few 
per cent. This "Pandemonium effect," originally pointed 
out in Ref. (3^ , was raised again recently (ssj with regard 
to superallowed /3 decay in p, /-shell nuclei. Following the 
approach advocated in these references, we use a shell- 
model calculation to compute the weak /? branches and 
include its predicted strength in the analysis. The model 
space used is the full s, d shell with the effective inter- 
actions USD of Wildenthal [sj] and the two more recent 
updates USD-A and USD-B of Brown and Richter [ist . 

We include in our analysis of the branches and yields a 
total of 51 excited states in '^^S. Our shell- model calcula- 
tion correctly predicts all of the 9 lowest (0, 1, 2)"*" states 
with < 7.2 MeV reported in Detraz et al. [2^. We 
find that the RMS deviations of the shell-model calcula- 
tion from the known excitation energies [l2j are quite 
good: 120 kcV (USD), 209 keV (USD-A) and 172 keV 
(USD-B). This is a gratifying indication that the shell 
model is performing well in this s, d-shell nucleus. Even 
though selection rules prohibit (3 decays to the six lowest 
(3,4)^ states, they are included in the analysis when ac- 
counting for 7-ray de-excitations. The shell- model calcu- 
lations identify approximately 40 (3 transitions to states 
whose excitation energies in '^^S lies between 7.485 and 
w 11.8 MeV. Unfortunately, the high density of states 
in this energy range makes a state-by-state comparison 
difficult, especially for the 2+ states. Based on the good 
correspondence of excitation energies and 7 de-excitation 
branches, we are able to identify 6 of the shell-model 
states in this region with ones in the ENSDF Data Ta- 
bles. None of the other 30 shell-model states individu- 
ally has a /3-transition strength greater than 0.3%, but 
cumulatively they sum to 0.50% in the USD, 0.69% in the 
USD-A, and 0.55% in the USD-B calculations. We include 
these weak /3 strengths and de-excitation 7 rays predicted 
by the shell model in our overall analysis. 

In the analysis, a /3 branch could be deduced as long as 
there is at least one decay 7 ray lying within the 7.35 MeV 
energy range of our HPGe. The ground-state branch and 
higher excitation-energy shell-model-state branches that 
were not observed in this experiment were included in 
the analysis as missing strength. For the ground state, 
we take the branch to be (1.0^0 5)% as determined by 
Armini et al. (25j , and the combination of all the unseen 
shell-model states at energies above 7.2 MeV is taken 
to be the average of the USD, USD-A and USD-B cal- 
culations with an uncertainty that spans the variation: 
(0.60±0.10)%. The a-particle and proton-emitting states 
in the 8.7 to 11.1 MeV excitation energy range reported 
by Honkanen et al. 3l| are not separately included be- 
cause their summed (3 strength of (0.080 ± 0.005)% is 
significantly less than and no doubt already included in 
the missing strength predicted by the shell model. 
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RESULTS 
Experimental Results 

The excitation energies in ^^S, the /3-decay branches 
and the log ft values determined in this work are shown 
in Tabic HI The states up to 7.2 MeV each have multi- 
ple 7 rays which were observed in this work. Thus the 
7 branches in these cases could be fit using the proce- 
dure described earlier. Furthermore, the fitting routine 
allowed us to treat the excitation energies of these states 
as free parameters, and fit the observed lines such that 
the Ex values minimized the of the calculated 7 ener- 
gies. For the highest three energy levels listed in Table HI 
only one de-excitation 7 ray was observed; therefore, in 
these cases the 7 branches had to be taken from previous 
work 12| , and the excitation energy had to be calculated 
solely on the one calibrated 7 ray that was observed. 

In Table HI] we list the 7 branches for the states which 
have multiple 7 rays within our 7.35 McV energy range. 
The 7 energy in this table is calculated by first averaging 
the ENSDF -Ex's with our own (the first two columns 
of TableHI and then calculating the Dopplcr-corrcctcd 
based on the new excitation energies. 

Figure [9] graphically depicts the results listed in Ta- 
bles [J and [ill The figure also includes within it the refer- 
ences to previous results which were used in the analysis. 
For example, the 7 ray from the 1'^ state at 9206 keV 
was outside the energy range of our HPGc, so we took 
its 7 de-excitation branching ratio from the ENSDF Data 
Tables in order to fit the /3 branch based on the observed 
7 ray to the 2^ first excited state. 

In both Tables |T] and |lTl the first uncertainty in the 
branches is statistical and the second is systematic. The 
sources of systematic error considered include: the cuts 
made on the data (the /3/HI ratio, the /3-7 timing win- 
dows for real versus accidental coincidences, and the 
start/stop of the counting time following transfer of the 
activity); the uncertainty in the 7 photopeak areas and 
total efficiencies; the (3 efficiency; the effective interac- 
tion of the shell model used for weakly-fed states above 
7.2 MeV; including/excluding the ENSDF 7 branches for 
higher levels (and varying them by their uncertainties 
when included); and the lo;5% uncertainty in the ground 
state branch. For the 7 branches in Table |IT1 the sum 
of the quoted probabilities for decay from a given ex- 
cited state in '^^S is not necessarily 100%, i.e. contrary to 
Eq. ([H, it may seem that J2j<iJi,j ^ 1; this is because 
of possible-but not statistically significant-peaks where 
the data only allows us to place limits on the branch. As 
with the 7 yields discussed below, a branch is quoted only 
if the area of a 7-ray peak in the HPGe was larger than 
its total uncertainty. If this condition is not satisfied, we 
instead quote an upper limit on the branch at the 90% 
C.L. For example, the 2^ state at 7115 keV has four 



statistically significant branches which have a combined 
probability of only 97.1 ± 5.3ti;2%. For the other two 
branches, the fit converged to results that were consis- 
tent with zero and for which only upper limits may be 
quoted: 1.2±2.3t° 5% (2+) and 1.7±2.0t^ j% (0+). All 
together, including these statistically insignificant transi- 
tion strengths, the total probability is 100%; the "miss- 
ing" 2.9% in Table M for the decay from the 2^ state 
is potentially within transitions which are below our de- 
tection sensitivity. Note, however, that the 1-a upper 
limit on the sum of the statistically significant branches 
includes 100%. 

With the /3 and 7 branches established, we are able to 
calculate the 7 yields. The reason we present the yields 
after the discussion of the branches is because of the small 
summing corrections which depend on the branching ra- 
tios. The results are listed in Table IIIII where we again 
compare the excitation energies to ENSDF [13] and the 
yields to the work of Detraz et al. [2^ . 

There is generally good agreement with the results of 
Detraz et al., although we find significantly less strength 
in the /3 branches to the 3.78 and 4.28 MeV levels. 
We attribute this discrepancy to the fact that many of 
the higher levels not considered in Ref. [2^ 7 de-excite 
through these levels; thus though our 7-ray yields for 
these states are in good agreement, our 7-ray feeding 
from higher levels results in a smaller deduced /3 branch 
for these states. Another difference from Detraz et al. is 
seen with the 7190 keV 1^ level where we see less than 
half as much 7 yield, and find a /3 branch that is 30% 
smaller. It is difficult to comment on this discrepancy 
since an efficiency curve is not provided in Ref. [29[. 

The shell-model predictions of the /? branches for 
the five highest energy levels of Table [1] are: 0.22(4)% 
(7536 keV); 0.05(4)% (7637 keV); 0.10(3)% (7921 keV); 
0.06(1)% (8125 keV); and 0.06(1)% (9208 keV). The 
agreement is quite reasonable, where the only signifi- 
cant difference (> 2a) seen is in the branch to the 8125- 
keV level; here the shell model calculation predicts w 3 x 
more strength than our limits allow for the branch to the 
8125-keV state. Note that the sum of these branches 
in the shell model, 0.49(7)%, is in perfect agreement 
with the corresponding sum of the observed branches: 
(0.51^0;^°)%. Given the high density of states in this 
energy range, it is a testament to the quality of the shell- 
model in this case that these branches are reproduced 
so well. It further justifies our use of the shell model to 
account for the Pandemonium effect as discussed earlier. 

In addition to reducing the uncertainties in all of the 
branches/yields by approximately an order of magnitude, 
we have observed 22 more 7 branches and three new (3 
branches compared to Detraz et al. [2^. For the ten 7 
and two (3 transitions where wc did not observe a statis- 
tically significant branch, a 90% CL limit is quoted. 
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TABLE I. Excitation energies in S, deduced /3 branches and the comparative half-lives of the decay of CI. The spin 
assignment of the 7637- and 7921-keV levels are determined from a comparison with our shell-model calculation. Not listed 
are the weakly-fed states between 7.4 — 11.8 MeV predicted by the shell-model to have a total branch of (0.60 ± 0.10)%. 
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TABLE II. 7-decay branches, in percent, for excited states in ^'^S and comparison to the currently accepted values in ENSDF [l^ . 
The 2230.48(14) keV 2^ transition to the Ot ground state is assumed to be 100% since there are no other known levels to 
which the 2^ state may decay. The energies of the 7 rays are calculated based on the adopted energy levels, i.e. the average 
of our work and ENSDF, the third column in Table H] and shown in Fig. m 
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FIG. 9. 13+ decay scheme for ^^Cl, summarizing the 7 and /3 branches in Tables |l] and |ll] determined from this work. All 
branches are expressed in percent. The excitation energies are the weighted average of the present work and the accepted 
values from ENSDF [IJ] (third column of Table 
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TABLE III. Yields of 7-rays following the /? decay of '^^Cl. All values are in percent, and limits correspond to the 90% 
confidence level. The overall normalization is fixed such that the transition to the ground state is consistent with l.Olg 5% as 
reported by Armini, et al. [23] and includes the (0.60 ± 0.10)% unobserved, weakly-fed shell-model states in the energy range 
of 7.4 - 11.8 MeV. 
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It 


(4695) 




0.027±0.008 


503.4 


±1.0 




2t 


(4282) - 


->ot 


(3778) 




< 0.04 


413.5 


±0.5 




It 


(4695) ' 


2+ 


(4282) 




< 0.019 



"Calculated from the adopted levels of Ref. [12] . 
''Ref. [29|. 



Comparison to Shell-Model Calculations 

Shell-model calculations have been performed for the 
states involved in the /3-decay of ^^Cl, which has a spin- 



parity of 1+. For transitions to 0+ and 2+ states in 
■^^S, the strength is pure Gamow- Teller whereas for the 
isobaric analogue transition, the decay is almost pure 
Fermi. Decays to non-analogue 1+ states can also in- 
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elude a Fermi component (via isospin mixing) and so for 
an experimental branch to one of these states, we can 
proceed in one of two ways: 

1. make some assumptions about the Fermi contribu- 
tion and deduce |M(GT)|. 

2. make some assumptions about the Gamow- Teller 
contribution and deduce |M(F)|. 

In the next section we discuss Gamow- Teller matrix 
elements after making some minimal assumptions about 
the Fermi contribution. We will compare our experimen- 
tal M(GT) values with shell-model computations to say 
something about the quality of the USD wave functions. 
Following this, we take the other approach by assuming 
that the USD values for B{GT) arc correct allowing us to 
deduce B(F). In turn we can calculate an experimental 
value for the isospin-mixing parameter, Sci, which we 
compare with theoretical calculations. 



Gamow-Teller matrix elements 

Here we present a comparison of calculated versus ex- 
perimental i?(GT)'s. Because several of the transitions 
are quite retarded, with log ft values exceeding 5 (see Ta- 
ble HI , the spectrum shape may depart significantly from 
the allowed shape. To proceed, we use a shell-model cal- 
culation to compute the shape correction function C(W) 
as described in the appendix of Ref. 3. We define an 
"exact" statistical rate function as 



Wo 



pW{Wo ~ Wf F{Z, W) C{W) dW, (6) 



where W = Ee/me is the electron total energy in elec- 
tron rest-mass units, Wq is the maximum value of W, 
p = {W^ — ly^^ is the electron momentum, Z is the 
charge of the daughter nucleus, and F{Z, W) is the Fermi 
function. The usual statistical rate function, /, as used 
for example in Table U to obtain log ft values, puts the 
shape correction function C{W) to unity. Taking the 
Fermi strength to be zero for these non-analogue tran- 
sitions, we obtain an experimental value for the GT 
strength from 



^cxp(GT) — 



2Tt 



/oxacti(l + S'ff) ' 



(7) 



and an associated experimental GT matrix element, 
Mexp(GT), defined by: 



Sexp(GT) = g2_,^|Me.p(GT)| 



(8) 



Because /exact depends on the GT matrix element via 
C(W), we use an iterative procedure which is explained 



below. In Eq. ^ we use J^t°*-^°'^ = 3071.81(83) s 
from the survey of Hardy and Towner [l|. The par- 
tial half- life, t, is calculated from the ^^Cl half- life of 
ti/2 — 298(1) ms [m and the branches, R, listed in Ta- 
blelH corrected for a small electron-capture fraction. Pec- 



.---fa 



Pt 



ECj 



(9) 



The correction S'^ is the transition-dependent part of 
the radiative correction and is obtained from a stan- 
dard QED calculation that depends on Z and W. It 
is evaluated to order a and Za^, with the order Z'^a^ 
terms estimated 



- 39[ . Thus with the quantities on the 
right-hand side of Eq. ([7| determined, an experimental 
-Soxp(GT) value is obtained. This relates to the Gamow- 
Teller matrix element, Mcxp(GT), as shown in Eq. ([S]). 
It requires knowledge of the ratio of the axial- vector to 
vector coupling constants, denoted gA, for which an ef- 
fective value is used, gA,cS in the context of shell-model 
calculations in finite model spaces. In the s. d shell, the 
systematic studies of Wildenthal and Brown 4^ 41 1 have 
shown the effective value to be of order unity, so we take 
gA.cS = 1- 

In the shell model, the calculation of the Gamow-Teller 
matrix element is based on: 



Af(GT)-^(/|ai&^K)(a|GT|, 



(10) 



a, 13 



where aj, creates a neutron in quantum state a, bp anni- 
hilates a proton in quantum state /3 and (q;|GT|/3) is the 
single-particle Gamow-Teller matrix element. The ma- 
trix element of a^^bp in the initial and final many-body 
states are known as the one-body density matrix elements 
(OBDME). The same OBDME used in the construction 
of the shape-correction function, C{W), are also used in 
the shell-model evaluation of M(GT). Our procedure 
is to scale one of these OBDME, recompute C{W) and 
/exact; and obtain a new Mexp(GT) from Eq. ([7]). We 
then repeat the procedure, refining the scaling at each 
step until the theory input matches the experimental out- 
put. We found that it did not matter which USD calcu- 
lation we started from, or which OBDME we scaled; the 
convergent result for the Gamow-Teller matrix element 
was always the same. Thus the method is quite stable. 

In Table IIVI we list the partial half-lives, t, the sta- 
tistical rate function /exact at convergence, the radiative 
correction (5^, and the deduced experimental Gamow- 
Teller matrix clement. We also list the theoretical values 
from three shell-model calculations with effective inter- 
actions USD, USD-A and USD-B, now without any of the 
adjustments to the OBDME discussed above. Thus, Ta- 
ble [IV] compares Mexp(GT) values with theoretical ex- 
pectations. The comparison is very favourable with the 
RMS difference w 0.07. Where the matrix element is 
quite large, say il/(GT) > 0.2, theory does exceedingly 
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TABLE IV. Deduced experimental Gamow- Teller matrix elements, Moxp(GT), in the decay of ^^Cl, and comparison with three 
shell-model calculations. The RMS deviations in the predicted |M(GT)| compared to experiment are 0.073 (USD), 0.068 (USD-A) 
and 0.067 (USD-B). 



Experiment 



Theory |A/(GT)| 



in ^'^S" 




T 


t[s] 


r b 
J exact 


3'r [%] 


|Af(GT)| 


USD 


USD-A 


USD-B 


9207.1 ±0.7 


1 + 


1 




133.2±0.2 


1.70(3) 


0.581^5 


0.300 


0.348 


0.298 


8125.4 ±0.2 




1 


> 150 


689.2±2.7 


1.55(3) 


< 0.20 


0.123 


0.146 


0.134 


7921.8 ±0.9 









879.5±1.3 


1.52(3) 


0.087±0.019 


0.130 


0.148 


0.183 


7637.0 ±1.0 







> 1200 


1261±4 


1.49(3) 


< 0.05 


0.121 


0.022 


0.056 


7535.6 ±0.9 




1 


161±18 


1392±3 


1.48(3) 


0.164±0.009 


0.193 


0.165 


0.192 


7190.3 ±1.1 


1^ 





48±4 


2024±4 


1.44(3) 


0.250±0.010'= 


0.143 


0.117 


0.162 


7115.1 ±0.8 


2^ 


1 


48±3 


2232.4±2.4 


1.44(3) 


0.238±0.008 


0.259 


0.166 


0.228 


7001.4 ±0.4 




1 


1.327±0.012 


2413. 0±1. 7 


1.42(3) 




0.012 


0.064 


0.036 


6665.9 ±0.8 


2+ 





14.3±0.5 


3411±4 


1.40(3) 


0.353±0.006 


0.428 


0.273 


0.365 


5548.4 ±0.8 







7.79±0.23 


8583±7 


1.30(3) 


0.301±0.004 


0.227 


0.333 


0.302 


4695.3 ±0.4 


It 





4.89±0.07 


15702±30 


1.24(3) 


0.281±0.002'= 


0.280 


0.309 


0.346 


4281.81±0.28 







13.7±0.5 


21280±10 


1.22(3) 


0.145±0.003 


0.071 


0.059 


0.085 


3778.3 ±0.7 







31±3 


28009±26 


1.18(3) 


0.083±0.003 


0.076 


0.137 


0.077 


2230.56±0.14 







0.504±0.003 


67470±40 


1.11(3) 


0.423±0.001 


0.423 


0.421 


0.457 


ground state 







30lf 


167900±1600 


0.94(3) 




0.004 


0.089 


0.024 



"The average of currently accepted values and the present work, i.e. the third column of Table U] 

''The error bar reflects both the uncertainty in the Q- value and the range obtained for difi'erent shell-model calculations 
of the shape-correction function, C{W). 

'^Foi these non-analogue 1^ transitions the experimental derived value is ^yW(Wy+^{GT) , where the Fermi contribution 
originates from isospin mixing as discussed in the text. For the 4695-keV level, the Fermi contribution is theoretically 
expected to be very small, so the value quoted is likely a good estimate for |Moxp(GT)|. For the 7190-keV level, the Fermi 
contribution is likely substantial, so a conservative estimate for |Moxp(GT)| is in the range < |Afoxp(GT)| < 0.260. 



well. For retarded transitions with M(GT) < 0.2, theory 
does not perform as well, but here the small values are 
a consequence of cancellations among shell-model ampli- 
tudes which are much harder to get precisely right. The 
same comparison is presented in Fig. [TOl where the inte- 
grated B(GT) values are displayed as a function of exci- 
tation eriergy. As has been found by Brown and Wilden- 
thal (40l.l4lj the USD effective interaction in s, d-shell nu- 
clei gives a reasonably accurate picture of Gamow- Teller 
properties in these nuclei. The newer interactions, USD-A 
and USD-B, perform equally well. 



Isospin-symmetry breaking in Fermi transitions 

We switch our attention to the /3 transition to the 7001- 
keV, lj,r=l isobaric analogue state (IAS). This tran- 
sition is a mix of Fermi and Gamow- Teller components; 
therefore from the partial half-life alone it is not possi- 
ble to deduce the Gamow- Teller matrix element. Hence 
the gap in Table |TVl However, the shell-model calcula- 
tion for the Gamow- Teller matrix element predicts for the 
IAS state a very small value indeed. This is a fortunate 
happenstance: it gives us the opportunity to study this 
transition as if it were a pure Fermi type, compare it with 
the precisely measured pure Fermi transitions between 
0"*" states, and deduce the amount of isospin-symmetry 




Excitation energy [MeV] 



FIG. 10. (Colour onhne) The observed integrated B{GT) 
as a function of excitation energy for the decay of ^^Cl and 
comparison to shell-model calculations. For the experimen- 
tal histogram (thick solid line with the shaded region rep- 
resenting its uncertainty bars), we have assumed zero Fermi 
strength in the 4695-keV state, zero Gamow- Teller strength 
in the IAS state, the range of Gamow- Teller strength given 
in the footnote of Table IIVI for the 7190-keV state, and in- 
cluded the /3-delayed particle emitting states from Ref. |31{ . 
The hatched region indicates strength that may be missed 
in the experiment because the phase space for P decay is too 
small resulting in /3-delayed branches less than 0.0015% which 
would not have been reported in Ref. [3l|]. 
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breaking (ISB) in this transition. A fairly large ISB effect 
is anticipated because in '^^S the IAS state is only 189 keV 
away from the 7190-keV state; a state with the same 
spin but different isospin. Perturbation theory predicts 
that when two states of the same spin are close together 
in the spectrum, Coulomb and charge-dependent nuclear 
forces induce a degree of isospin-symmetry breaking that 
is inversely proportional to the square of the energy sep- 
aration of the two states. For a separation of 189 keV, 
mixing at the several percent level can be anticipated. 

The partial half-life for decay to the IAS state has been 
determined with 1% accuracy, namely t^^^ = 1.331(12) s. 
From this, the ISB correction 6c can be determined to 
15% accuracy from the equation: 

U.ctt'^'^il + 6'ij){l + S^s - 6c) 

K 

~ Gl{l + ^l)[B{V) + B{GT)] (11) 

~ B(F) + S(GT)' 

Here Kj(h,cf = 27r^?iln2/(meC^)^ is a constant and 
Gy is the vector coupling constant characterizing the 
strength of the vector weak interaction. The quan- 
tity KIG\(\ -h A^) is taken from the precision work 
on 0"*"— > 0"*" superallowed transitions and is expressed in 
terms of Tt^ introduced in Eq. ([7]). The radiative 
correction has been split into three pieces: (a) a nucleus- 
independent term, A]^, is included in Ti^^'^^ ; (b) a 
trivially nucleus-dependent term, (5^, is calculated to be 
1.421(32)%; and (c) a second nucleus-dependent term, 
(5ns 1 is small but requires a nuclear-structure calculation 
to be evaluated. It is convenient to place (5ns S'lid be 
together as both are dependent on shell-model nuclear- 
structure calculations. Finally i?(F) is the square of the 
Fermi matrix element, i3(F) = |A/oP = 2 for T=l tran- 
sitions in the isospin-symmetry limit, and i3(GT) is the 
square of the Gamow- Teller matrix element, Eq. ([7]). For 
i3(GT), we take the three theoretical values from the 
shell-model calculation using USD, USD-A and USD-B ef- 
fective interactions, average them and assign an uncer- 
tainty equal to half the spread between the largest and 
smallest calculated values: S(GT) = 0.002 ± 0.002. On 
rearranging Eq. (|11[) . we obtain 

_e , , 2(J-t°"-^"") 

'^^s^-P-^ /exact tiAS(l + <5y[i3(F) + S(GT)] 
= 5.4(9)%, (12) 

a substantial isospin symmetry breaking term, the largest 
yet determined in a superallowed Fermi transition. 

In what follows we present a shell-model calculation of 
6c — ^NS following the procedures developed by Towner 
and Hardy [3|. First, for the ISB correction 6c defined in 
Eq. ([1]) , the technique is to introduce Coulomb and other 



charge-dependent terms into the shell-model Hamilto- 
nian. However, because the Coulomb force is long range, 
the shell-model space has to be very large indeed to in- 
clude all the potential states with which the Coulomb 
interaction might connect. Currently this is not a practi- 
cal proposition. To proceed, Towner and Hardy 01 divide 
(5(7 into two parts: 

6c — 6c\ + 6c2- (13) 

For (5(71, we perform a shell-model calculation in the trun- 
cated ^hbj model space of the s, (i-shell orbitals. Charge- 
dependent terms are added to the charge-independent 
Hamiltonians of USD, USD-A and USD-B. The strengths of 
these charge-dependent terms are adjusted to reproduce 
the h = -5.4872(35) MeV and c = 0.1953(37) MeV [13 
coefficients of the isobaric multiplet mass equation 
(IMME) as applied to the 1+, r= 1 states in A = 32, the 
triplet of states involved in the /3-transition under study. 
As mentioned already, the bulk of the isospin mixing in 
the IAS occurs with the neighbouring 1;^ state. This ob- 
servation is used to constrain and refine the calculation. 
In the limit of two-state mixing, perturbation theory in- 
dicates that 

6cx oc 1/(A£;)2, (14) 

where Ai? is the energy separation of the analogue and 
non-analogue 1"*" states. Thus it is important that the 
shell-model Hamiltonian produce a good-quality spec- 
trum of 1+ states. The shell model calculation has vary- 
ing degrees of success in this regard. For the 1+ states 
in A = 32, the separation between the IAS and the third 
i;[, r=:0 state is observed to be 188.9±1.2 keV. The shell 
model calculates this separation to be 184 keV with USD, 
248 keV with USD-A and 387 keV with USD-B interactions. 
These are quite respectable results given the inherent ac- 
curacy of a shell-model calculation for predicting ener- 
gies. However, for a reliable 6c\ calculation, this spread 
in Ai? values is quite a problem. To cope with this, the 
Towner-Hardy recommended procedure is to scale the 
calculated 6c\ value by a factor of (Ai?)jjjpQ/(Ai?)^xp) 
the ratio of the square of the energy separation of the 1"*" 
states in the model calculation to that known experimen- 
tally. After this is done, the 6c\ values obtained in the 
three shell-model calculations are reasonably consistent: 
6ci = 3.73% for USD, 3.32% for USD-A, and 4.19% for 
USD-B. We average these three results and assign an un- 
certainty equal to half the spread between them to arrive 
at: 

6cx = 3.75(45)%. (15) 

For the calculation of 6ci we need to consider mix- 
ing with states outside the OSw shell-model space. The 
principal mixing is with states that have one more radial 
node. Such mixing effectively changes the radial func- 
tion of the proton involved in the /3 decay relative to 
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that of the neutron. The practical calculation, therefore, 
involves computing radial overlap integrals with modeled 
proton and neutron radial functions. Details of how this 
is done are given in Ref. [^. The radial functions arc 
taken to be eigenfunctions of a Saxon- Woods potential 
whose strength is adjusted so that the asymptotic form 
of the radial function has the correct dependence on the 
separation energy. The initial and final A-body states are 
expanded in a complete set of (A— l)-parent states. The 
separation energies are the energy differences between 
the A-body state and the (A— l)-body parent states. A 
shell-model calculation is required to give the spectrum 
of parent states and the spectroscopic amplitudes of the 
expansion. For the three USD interactions, we compute 
5c2 = 0.827% for USD and 0.865% for both USD-A and 
USD-B. Our adopted value is: 
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FIG. 11. Our determination of the isospin-symmetry- 
breaking correction for ^^Cl (filled circle) and calculations for 
^^Cl as well as other superallowed transitions (open points). 
The three components, 5ci, 5c2 and 5ns, are shown sepa- 
rately. The measurement and prediction for ^^Cl, particularly 
the Sci component, is significantly larger than in any of the 
0^^- 0^ transitions. 



The uncertainty, calculated in the same manner as de- 
scribed in Ref. (j) , represents the range of results for the 
USD interactions, the different methodologies considered 
in adjusting the strength of the Saxon- Woods potential, 
and the uncertainty in the Saxon- Woods radius parame- 
ter which was fitted to the experimental charge radius of 

Finally, we need an evaluation of the nuclear-structurc- 
dependent piece of the radiative correction, (5ns- Such 
a term arises because in a many-body system such as 
a nucleus, the electromagnetic interaction and the weak 
interaction that collectively induce a radiative correction 
do not have to interact with the same nucleon in the nu- 
cleus. When these interactions occur with different nu- 
cleons, the process is described by two-body operators. 
The evaluation of matrix elements of two-body operators 
depends in detail on the nuclear structure of the states in- 
volved. Such calculations were first made in 1992 43l l44| 
and updated two years later [i^. We follow the latter 
reference and compute (5ns for each of the s, ci-shell ef- 
fective USD interactions. Essentially the same result was 
obtained in each case. We adopt the average value of 



3ns 



-0.15(2) 



(17) 



The result is a very small correction, about 3 times 
smaller than the uncertainty in 5c- 

Adding together Eqs. UHl), (US]) and (H?]), we obtain 



5c - <5ns = 4.8(5)%, 



(18) 



which agrees with the experimental result of 5.4(9)% 
of Eq. (|12p within stated uncertainties. A comparison 



of this result with calculations for the 0+— 0"^ cases is 
shown in Fig. [TTJ As one can clearly see, the correction 
in 32 CI is about five times larger than the typical < 1% 
values found for the s,c?-shell nuclei in 0^— > 0^ super- 
allowed transitions. The TH model, which has already 



been shown [46| to reproduce the nucleus-to-nucleus vari- 
ation of ISB effects in superallowed zerotozero transitions 
required by the CVC hypothesis, is shown here to pro- 
duce a much larger ISB effect as again verified by the 
current experiment. 

Let us now briefiy consider the isospin mixing in the 
non-analogue 1+ states at 7190 and 4695 keV, and de- 
duce experimental values for the degree of mixing in the 
same way as just described for the IAS state. We assume 
the shell-model calculation for B{GT) to be correct, as- 
signing an uncertainty equal to half the spread between 
the different results obtained with the USD, USD-A and 
USD-B interactions. For the 7190-keV state, the shell- 
model calculations yield S(GT) = 0.020(6), and for the 
4695-keV state B(GT) = 0.098(21). An experimental 
value for B{¥) -\- _B(GT) is computed from a rearranged 
Eq. (Hg): 



Scxp(F) 4- Bcxp(GT) 



2J-tO 



/cxacti(l + l5p) ' 



(19) 



where the nuclear-structure-dcpendcnt radiative correc- 
tion, (5ns, is ignored. Inserting the experimental values 
into the right-hand side of Eq. ([TO]) , we obtain 



5cxp(F) + Sexp(GT) = 0.062 ± 0.005 
-Bcxp(F) = 0.042 ±0.008 

for the 7190-keV state, and 



(20) 



Bexp(F) + Sexp(GT) = 0.0790 ± 0.0012 
-Bcxp(F) = -0.019 ±0.021 
or Sexp(F) < 0.014 (21) 



for the 4695-keV state. In the limit of exact isospin 
symmetry the B{Y) values would be zero for these non- 
analogue transitions. So the non-zero value in Eq. (|20p 



16 



is a further indication that isospin-symmetry breaking is 
present. For non-analogue transitions, we define 



B{F) 



|A/oP'5£i(l 



(22) 



where S^ji is the isospin-symmetry breaking correction 
for the n*'^ 1+ state in '^^S computed in the s, d shell- 
model space, and dc2 the radial overlap correction rep- 
resenting Coulomb mixing beyond the Ohuj model space. 
In Eq. (|22p . we drop the 6ci5c'2 cross term as being neg- 
ligible. From the experimental results in Eqs. ([20|l and 
([2T|l . we determine (^^i^^xp = 2.1(4)% for the 7190-keV 
state and (5^i_exp < 0-7% for the 4695-keV state. The 
theory calculation that produced the result for dci in 
Eq. (fT5|) also gives as a by-product values of (5^-^ for the 
non-analogue transitions. For the 7190-keV and 4695- 
keV states, theory predicts 6^^ — 3.2(3)% and 6^^ = 
0.04(1)% respectively. Evidently there is a small discrep- 
ancy between theory and experiment for the symmetry- 
breaking in the 7190-keV transition. On the theory side, 
if we accept the symmetry-breaking in the IAS state to 
be correct, then it is most likely to be correct in the 
7190-keV state as well. This is because the situation is 
close to 2-state mixing, and the loss of Fermi strength to 
the IAS state is recovered in the 7190-keV state. On the 
experimental side, the result depends on the correctness 
of the Gamow- Teller strength calculated with USD wave 
functions. If the B{GT) were over-estimated by the USD 
calculation, then theoretical and experimental values for 
the isospin symmetry mixing in the 7190-keV state could 
easily be reconciled. 



CONCLUSIONS 

We have measured relative 7-ray intensities and de- 
duced /3-decay branches for the decay of ^^Cl. We have 
observed 3 new P branches, 22 new 7 lines, placed limits 
on 2 other /3 branches and 10 other 7 transitions, and 
have improved the precision on previously known yields 
and branches by about an order of magnitude. 

In total, twelve /3 branches have been measured in the 
decay of '^^Cl. Eleven of these are Gamow- Teller transi- 
tions and one is predominantly Fermi. For the Gamow- 
Teller transitions, the GT matrix element has been de- 
termined and compares favourably with shell-model cal- 
culations using USD effective interactions. These calcu- 
lations also find the Gamow- Teller component in the 
IAS transition to be very small, indicating this tran- 
sition is almost pure Fermi-like. Thus, this transition 
can be analyzed in an identical way to that used for the 
0+— > 0+ superallowed transitions. We extract a sizable 
isospin symmetry breaking correction for this transition, 
Sc — ^NS = 5.4(9)%, which agrees well with a theoretical 
value of 4.8(5)%. 

In addition, the improved precision in the relative 7- 
ray intensities can be used for a more precise determina- 



tion of 7-ray efficiencies in the decay of ^^Ar [l3|. The 7 
intensity from the lowest T = 2 state in '^^Cl is of interest 
for measuring isospin symmetry breaking in the T ~ 2 
superallowed decay of '^^Ar. Presently the 7-decay in- 
tensities from the decay of '^^Ar are limited by statistical 
precision, but the present work opens the possibility of 
determining its 7 branches to higher precision in future 
experiments. 
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